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Abstract 

In this paper we compute the Casimir energy for a coupled fermion-pseudoscalar field system. 
In the model considered in this paper the pseudoscalar field is static and prescribed with two 
C*) ' adjustable parameters. These parameters determine the values of the field at infinity (±#o) and 

. its scale of variation (/x) . One can build up a field configuration with arbitrary topological charge 

■ by changing 9o, and interpolate between the extreme adiabatic and non-adiabatic regimes by 

changing /x. This system is exactly solvable and therefore we compute the Casimir energy exactly 
and unambiguously by using an energy density subtraction scheme. We show that in general the 
Casimir energy goes to zero in the extreme adiabatic limit, and in the extreme non-adiabatic limit 
when the asymptotic values of the pseudoscalar field properly correspond to a configuration with 
an arbitrary topological charge. Moreover, in general the Casimir energy is always positive and on 
the average an increasing function of 8q and always has local maxima when there is a zero mode, 
showing that these configurations are energetically unfavorable. We also compute and display the 
energy densities associated with the spectral deficiencies in both of the continua, and those of the 
bound states. We show that the energy densities associated with the distortion of the spectrum of 
the states with E > and E < are mirror images of each other. We also compute and display 
the Casimir energy density. Finally we compute the energy of a system consisting of a soliton and 
a valance electron and show that the Casimir energy of the system is comparable with the binding 
energy. 

1 Introduction 

Dutch physicist Hendrick Brugt Gerhard Casimir published his famous paper [T] in 1948, where he 
found a simple yet profound explanation for the retarded van der Waals interaction between polarizable 
molecules. He showed that there exists a net force between grounded infinite metallic plates, placed 
a few micrometers apart, in a vacuum without any external electromagnetic field. The Casimir effect 
is essentially a direct consequence of the change in the vacuum energy of a quantum field when some 
boundary condition or background field is imposed. Since the Casimir's work, many papers have been 
written on the Casimir energy or the resulting force due to the presence of boundary conditions. They 
have calculated this effect for parallel plates [2], cylinders [3J, and spheres [I] (first considered by 
Casimir [5], DeRaad and Milton [5], and Boyer [7], respectively), and many other geometries [5]. In 
addition, many authors have considered this effect for various fields, different boundary conditions, 
dimensions and topologies [5]. There are various methods to calculate the Casimir energy, e.g. the 
Green function formalism [TD] , the heat- kernel method [TT] and the multiple-scattering expansions [12] . 
Also, as mentioned before, the presence of the nontrivial background fields such as solitons affect the 
zero-point energy. For example, to compute the quantum corrections to the mass of the soliton, the 
change in the quantum fluctuations of the vacuum energy caused by the presence of the soliton, should 
be taken into account. In 1974, Dashen et al [13] computed the one-loop correction to the mass of the 
kink of the 4 theory for the first time. After this work, several authors have used various methods to 
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compute the corrections to the soliton mass, such as the scattering phase shift method, zeta function 
analytic continuation technique, dimensional regularization technique, etc. [141 115] . Moreover, the 
Casimir effect appears in supersymmetric models. To investigate the validity of the BPS saturation by 
supersymmetric solitons, one should compute the change in the fluctuations of the vacuum energy to 
obtain the quantum corrections to the soliton mass (see for example 16 J. 

One of the first experimental attempts to verify the existence of the Casimir effect was conducted 
by Marcus Sparnaay for two parallel metallic plates in 1958 [T7]. The results was not in contradiction 
with the Casimir theory, but had a very poor accuracy. In 1997, using a plate and a spherical lens 
covered by a metallic sheet, Steve K. Lamoreaux measured the Casimir energy with high accuracy 
[18 . It was the first "successful" experiment to verify the Casimir effect. Since then, many different 
experiments have been performed to measure the casimir effects for various geometries, such as two 
parallel plates and a sphere in front of a plane, with ever increasing accuracy (see for example |19j). 

In this paper we consider a pseudoscalar field as a prescribed, spatially varying background field 
which is coupled chirally to a Fermi field. The background field in our model has two adjustable 
parameters, one related to its scale of variation and the second to its asymptotic value at spatial 
infinity. By varying the latter from zero to mr, we can build a field configuration with topological 
charge n. As it is well-known the nomenclature "soliton" refers to a static field configuration which is 
a solution to the classical field equations with finite energy and localized energy density. Topological 
solitons by definition carry nonzero topological charge. In a field theory in which this background field 
interacts with a fermion, such as the one considered in this paper, the exact shape of the soliton depends 
on the back-reaction of the Fermi field, which in turn depends on the particular fermionic state to which 
the soliton is coupled. Therefore, the soliton can in principle acquire infinitely many different shapes. 
However, since the back-reaction is usually small, all these different shapes are close to each other. 
However, we have obtained the exact (to within numerical approximation) non-perturbative solutions 
for the coupled system treating both fields as dynamical, verifying the perturbative result that the 
back- reaction is small [20 . The problem of obtaining the Casimir energy for any of these particular 
field configurations is also not exactly solvable. Therefore, in this paper we choose a simple field 
configuration which has proper topological charge and renders the problem exactly solvable. Although 
the field configuration that we choose is not a solution to the field equations, it possesses all other 
properties of a proper soliton. We shall henceforth refer to this configuration as a soliton for brevity 
of notation. We expect the exact Casimir energy for this topologically non-trivial configuration to 
illuminate many of the qualitative features of the problem with an exact soliton. Since this model is 
exactly solvable and we have the complete spectrum of the Fermi field [2T] , we are able to calculate the 
vacuum expectation value of the Hamiltonian directly and exactly for all values of the aforementioned 
parameters. We obtain the Casimir energy for this system by subtracting the zero-point energy in 
the presence and absence of the background field. In Section II we introduce the model and present 
a derivation for the expression of the Casimir energy. In Section III we obtain an exact expression 
for the Casimir energy and its density for our model. Then, we thoroughly explore the properties 
and behavior of the Casimir energy as a function of the parameters of our model. First, we depict 
the Casimir energy as a function of the scale of variation of the soliton with winding number one. In 
general there is a sharp maximum whenever the Fermi field has a zero-mode bound state. We obtain 
the expected results that the casimir energy tends to zero in the limits of extremely wide and sharp 
solitons. Second, we depict the Casimir energy as a function of Qq. In general -Ecasimir is an increasing 
function of 9q with very mild maxima in the form of cusps whenever there is an energy level crossing 
zero. Then, we exhibit the energy densities of the bound and continuum states and their sum for states 
with E < and E > 0, separately. These two energy densities are exact mirror images of each other, 
and from this we compute the Casimir energy density. In Section IV we summarize and discuss the 
results. 
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2 A fermion-pseudoscalar field model and the Casimir energy 



Consider the following Lagrangian describing a fermion interacting non-linearly with a pseudoscalar 
field 4> (x), in (1 + 1) dimensions, 



me 



(1) 



In order to have an exactly solvable model, Gousheh and Mobilia [21] chose </> (x) to be piecewise linear 
as follows 

{— #o for x ^ — I, 
ixx for — I ^ x ^ I, (2) 
+6<o for I sC x, 

where fj, = 9q/1 is the slope in the middle region. One could interpolate between the extreme adiabatic 
and non-adiabatic regimes by changing fi from zero to infinity, respectively. When 9q = mr, the soliton 
has winding number n. Figure ([T|) shows the form of <p (x). As mentioned earlier this simple topological 
configuration is used as an approximation to the kink of the Ac/) 4 theory. Although our simple static 
configuration is not a solution to the Euler-Lagrange equation, its total energy can be easily obtained 
from the corresponding Hamiltonian as follows 



U[(t>{x)]dx = 



m 



dx = ^|^o, 
15 



(3) 



where = 0q, the value of 4>{x) when x — > oo, and = /i, the slop' 



e at x = 0. 



<p(x) 



slope = jti 



Figure 1: The form of the background pseudoscalar field <j>(x). 



The procedure is to find the spectrum of the Fermi field in the presence of this background, for arbitrary 
values of /x and 9q. Then, for fixed /i, one could observe the changes in the spectrum of the fermion as 
a soliton with the topological charge n is built by changing 0q from zero to nix. This method was first 
introduced by MacKenzie and Wilczek [22] and is based upon the adiabatic method first introduced by 
Goldstone and Wilczek [23]. Since this model is exactly solvable, the complete spectrum of the Fermi 
field, including the eigenfunctions and their corresponding eigenenergies, has been obtained exactly 
for arbitrary values of fi and #o, and are given in |21j . Here we do not present the solutions in detail 
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Figure 2: The energies of the bound states as a function of 8q for jj, — 4m obtained from ™ sm(CO g o _ 
derived in [21] . Solid (dashed) lines show the bound states with positive (negative) parity. 



e = e -Wf) = (^+^\ (4) 



and the interested reader is referred to |21j . If we denote the Fermi field in (1 + 1) dimensions by 
ip = ( 1 ] , then we can define 

W2 

£ = e~ iEt ( , , 

and solve the equations of motion in terms of £. Since the Hamiltonian commutes with the parity 
operator, the eigenfunctions of the Hamiltonian can be chosen to be eigenfunctions of the parity 
operator as well. In the first representation (7 = <j\, 7 1 = 1173, 7 5 = 7°7 1 = 02)1 the parity operation 
is given by P4>{x : t) = a-y^) (— x, t). In the new representation i.e. Eq. (j4]) (7 = a%, 7 1 = ia±, 
j 5 — 7°7 1 = 173), it becomes P£ (x, t) — — <T2^ (— x, t). Figure (J2]) shows the energies of the bound 
states as a function of 8q for fj, = Am as an example. As shown in the figure, the bound states have 
alternating parity signs, starting with the positive sign. This graph will be important for the upcoming 
discussions. 

In order to obtain the Casimir energy, we should subtract the zero point energy in the absence 
from the presence of the background field. Therefore, the vacuum energies for both cases should be 
calculated. MacKenzie and Wilczek [35] assumed that the solutions in the presence of the soliton are 
complete. Later on Gousheh and Mobilia proved this assumption. Therefore, the Fermi field operator 
can be expanded in terms of either the modes in the absence or in the presence of the background field, 
as follows 



ip(x,t) = / — \b k u k (x,t) +dlv k {x,t) 
J-00 27r 1 J 



where Uk (vk) are the positive (negative) energy continuum eigenstates in the absence of the background 
field. Also, fj? (y^) are the positive (negative) energy continuum eigenstates and xib, (X26 t ) ar e the 
positive (negative) energy bound states in the presence of the soliton. Here, J2i refers to the possibility 
of multiple bound states and the superscript j refers to the parity of the continuum states. MacKenzie 
and Wilczek chose the coefficients of both positive and negative energy bound states to be annihilation 
operators. However, we choose annihilation (creation) operator for the positive (negative) energy 
bound states, as we do for the continuum states. We shall see one of the advantages of this choice 
when we compute the Casimir energy. In general the Hamiltonian can be expressed as follows 

/•+00 

H= dx [^(x,t)Hip(x,t)] . (6) 



Substituting the expression for the field operator tjj in Eq. (|5|) in the presence of the background field 
into Eq. ©, we obtain 



/+oo f r + OG J _ 

dx / / E [<^ f (M) + + E 4*1* (M) + /ixU*>*) 

/; 4 / EJ ^ xl „, ) ._,,; ( ,„ 1 

u n— ± 

[$L*a?iXUn *) + K'ounJh^bt (x, t)] }, (7) 



+ E 

i 

where, the ± superscripts on fvj, ound denote the sign of the bound state. Upon calculating the ex- 
pectation value of the Hamiltonian in the vacuum state, all the terms vanish except the terms with 
an annihilation operator on the left side and the corresponding creation operator on the right side. 
There are only two terms in this form which can be written in the following form by the use of the 
anticommutation relations 

4<t = - c fc> p + S ^S qp , /J/ = -flu + S u . (8) 
Therefore, we obtain the following expression for the vacuum expectation value of the Hamiltonian 

r>+oo />+oo 



dx I £ E (-Vp r +^) yp(x,t)4{x,t) 

/+oo 
^E Abound) xl bi (x,t) X 2 bi (x,t), (9) 
-°° i 

where | Q) is the vacuum state in the presence of the background field. From this expression we conclude 
that the entire contribution to the vacuum energy for fermionic systems originates from the negative 
energy eigenstates. Note that our separation of the bound states into the positive and negative energy 
ones has the advantage of the automatic inclusion of the negative energy bound states in the expression 
for the Casimir energy. To obtain the Casimir energy, one should subtract the zero point energy in the 
absence from the presence of the background field. Thus, we obtain 

f'Casimir — 

/ + oo /-+00 j n+OO 

_ dX l £ E (->/F+rt) 4H + J ^ (^und) X 2bi X2b t 

(-Vfc 2 +m 2 ) v\v k . (10) 



+°° , f +oc dk 

dx , 

2tt 



— oo 



Note that our problem possesses charge conjugation symmetry and therefore our expression for the 
Casimir energy given in Eq. (|10p is equivalent to the conventional one, where one would sum over all 
modes symmetrically with a factor of 1/2, while preserving the sign. 



3 The Casimir energy and its density for this model 

By substituting the expressions for the eigenstates into Eq. (llOj) . the Casimir energy can be calculated. 
Our prescription for subtracting two quadratically divergent integrals appearing in this equation, is to 
subtract the integrands with corresponding values of p and k, and then performing the integral. The 
change in the energy densities of the Dirac sea (E ^ — m) and sky (E ^ m) (i.e. the difference between 
the energy densities in the presence and absence of the background field) are 

£sk y, S ca(x) = / / ]T (E) ${X)$(X) ~ / - (E hcc ) g*{x)£» (*), (11) 

JO 2lT ~L J-oo 27T 
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where E = ± y/p 2 + to 2 , Ef YCC = ±\/k 2 + m 2 (± are for the Dirac sky and sea, respectively), = 
nUvi) for interacting Dirac sky (sea) and £| ree = Uk(vk) for free Dirac sky (sea). The bound state 
energy densities are 

£bound±(a;) = (Abound) ^W^Wi (12) 

i 

where the ± refers to the sign of the bound state energies, and £f, 4 = Xib;(X2t>i) for the positive 
(negative) bound states. Substituting the explicit expressions for the states given in [JTj into Eqs. ([TT]) 
and (p~2|> . we obtain the energy densities of the two continua and the bound states, shown bellow, 
for |x| I: 

£sky, S ca(x) = J -^EU\N c+ \ 2 + \N c _\ 2 )m[-cos[2p(\x\-l)](\ki\ 2 -p 2 ) 

x cosh [Im (C) I] + sin [2p(\x\ - I)} 2plm(fci) sinh [Im (C) I] ] 

+ (\N c+ \ 2 - |iV c _| 2 ) [cos [2p(\x\ - I)} [E (|fc!| 2 +p 2 ) + 2p 2 Re(fci)] cos [Re (C) I] 

- sin[2p(|2;| -0]p[|fci| 2 +p 2 + 2^Re(fci)]sin[Re(C)/]]|, (13) 

£b ou„d ± (x) = 2AT 2 i? b ± ound e-^^l-'), (14) 

and for \x\ ^ 

£sky, S ea(a;) = jf - 1 + ( | ^V c+ 1 2 + |7V C _| 2 ) [l^ 2 + p 2 + 2ERc (ki) + 2m 2 ] 

x cosh [Im (0 X] +2m[£ + Re(fci)] cos [Re (C) x] (|^ c+ | 2 - |iV c _| 2 ) |, (15) 

6bo»nd±(x) = ^%^ { (Abound " 2A 2 ) COS [Re (C) (x I)} 

[Re (C)] I 

- Re (C) A sin [Re (C) {x - I)] + \y 2 + /xB^}, (16) 

where A = yjm 2 — E 2 ound , ( = ^//x 2 — 4 (A 2 — (where £7' = {£", inbound} depending on the case), 
fci : 2 = 5 (/i ± C) and the normalization constants are 



* = 4 



A [Re(0] 



Re(C) 

1/2 

^ (M + 2^^) 



. P± _ o \ 2 

bound ■ sin [2ZRe (£)] — A cos [2ZRe (£)] + A 



iV c± = \ cosh [Im (C) I 

*l| a _ 1 ^- 1/a 



A-i I 2 

11 1 ^ " ■ =p cos [Re (£) i] 2mE 



, ,1 2£^ + 4£Tte(fci) 



P 2 



All of the spatial integrals can be performed analytically. However, the remaining integrals over p 
in Eqs. (|13I15|) cannot be done analytically and we integrate them numerically. Figure Q shows the 
Casimir energy as a function of y, at 9q = it, that is a soliton of winding number one. 



This graph has a sharp maximum at the critical value /i rs 2.957to, where a bound state crosses E = 
precisely at 9q = tt. This crossing is similar to what is depicted in Fig. ([2]) for /i — Am where one bound 
state crosses E = at 6q ss 0.77-7T. The appearance of the maximum at the critical point implies that 
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Figure 3: The Casimir energy as a function of p, the scale of variation of the soliton, at #o = 7T. The points 
displayed are the results of the numerical integration and the solid line represents the best fit to the points. 
In the illustrated box we focus on small values of p to show the details of the maximum. This maximum is 
precisely a cusp. 



constructing a soliton with this slope needs more energy than the ones with other slopes. According to 
Jackiw and Rebbi's profound work [53], when the system has particle conjugation symmetry and there 
is a zero-energy fermionic mode, the soliton carries fermion number In our problem which lacks 
this symmetry, our analysis shows that this point is energetically unfavorable. However, this property 
might be shared by problems possessing this symmetry. 

One can easily obtain the behavior of the Casimir energy at 9q = 7r numerically for large values of /i 
and the result is " 0.00163908 + 18.5905/x - 52.6919/x 125 + 67.2916/x 1 5 - 47.1817/x 175 + 15.3126/r 2 
", where x = fi/m. Therefore, we can conclude that to within our numerical accuracy, it tends to zero 
as fx — > oo. This is the expected result, as we shall describe shortly. Gousheh and Mobilia obtained 
the difference between the total number of levels in the Dirac sea in the presence and absence of the 
soliton, and in the limit /i — > oo their result is as follows 

Q _ = _ 1 + r + °° dp m 2 sin(fl )cos(flo) 

2 Jo 7T _|_ m 2 [p2 _|_ m 2 sm 2 (# )] 

Using this relation, the Casimir energy for the case fi — > oo would be 

_m f + °° dp ro 2 sin(6> ) cos(6> ) ^ , , > 

2 JO n [P + m Sin (Po)\ i 

The last term appears automatically due to our choice of the expansion of the Fermi field in the 
presence of the background field (we choose the annihilation (creation) operator for positive (negative) 
bound states). When /i — > oo one bound state just reaches the Dirac sea at 0q = nn, i.e. it becomes 
a threshold bound state. Therefore, like all threshold bound states in one spatial dimension, it should 
be counted as one half bound state. Since its energy is — m, the last term in Eq. (I18[) is — m/2. 
Therefore, the Casimir energy is zero in this limit at #o = n7T i since the integrand term vanishes 
identically there. For small values of p. the graph of the Casimir energy can be approximated by 
"0.0000252628 + 0.252951x - 0.00272887a: 1 - 5 - 0.015381x 2 " . This shows that the Casimir energy tends 
to zero, to within our numerical precision as fi — > 0. Therefore, as the slope of the soliton decreases 
to zero, the vacuum energy approaches that of the trivial vacuum, in spite of the residual non-trivial 
boundary conditions. These two limits occur at any 8q — mr i.e. when we have proper solitons with 
winding number n. However, for arbitrary values of 9q, although the result -Ecasimir — > as fi — > still 
holds, for fi — > oo the Casimir energy is in general non-zero and has the following form 



-^Casimir ~ ^ 



1 + 




A" 


^ cos 6*o 






7T 



mcos6> e(-cos(9o)e (e - — A, (19) 
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where 9(x) denotes the usual step function defined by, 



9 (a;) = { for x < ; i for x = ; 1 for x > } , (20 s 



and [^-] extracts the integer part of -k. In general the jumps in the energy of continuum and bound 
states at the points where threshold bound states appear cancel out so that the Casimir energy is 
continuous. 

In Fig. (j3J) we show the Casimir energy as a function of 6q for /.i = 4m. Note that the Casimir 
energy is, on the average, an increasing function of 9q for fixed /i. Comparing Fig. (j4|) with Fig. ([2]), we 
can conclude that there are local maxima in the form of mild cusps in the graph of the Casimir energy 
when bound states cross E = 0, as we vary 6> . 

In Eq. (|10p the integrand of the spatial integral represents the Casimir energy density. In what 
follows we shall compute and display the changes for both the positive and negative parts of the 
spectrum for comparison. In general the change in the energy density naturally divides into two pieces. 
First is the contribution coming from the continuum states (\E\ to) and second the contribution 
coming from the bound states (\E\ < to). 




Figure 4: The Casimir energy as a function of #o for /i = 4m. The mild cusps in the graph are actually mild 
maxima and occur when bound states cross E = 0. This also occurs in Fig. ©, except that the maximum 
there is much more pronounced. 



The left graph in Fig. ([5]) shows the energy densities of negative and positive bound and continuum 
states for the case /i = 4to at 8q = ir, given in Eqs. (|13I14I15I16|) . Note that an overall positive energy 
density in the sea and overall negative energy density in the sky indicate that both of the continua 
contain spectral deficiency as compared to the free case. The sum of these deficiencies is precisely 
canceled by the sum of bound state densities, at each point of space [21 . The sum of e sca (x) and 
Abound- (x) (solid and dashed lines in the left graph, respectively) is the total distortion of the states 
with E < 0, and is depicted in the right graph in Fig. ([5]) by the solid line. Also, in this graph we show 
the sum of £ s k y (x) and Ebound+(^) (dotted and dotdashcd lines in the left graph, respectively), which 
is the total distortion of the states with E > 0, by a dashed line. This figure shows that these two 
total densities are exactly the mirror images of each other. It is important to note that in problems 
which possess charge conjugation symmetry, such as ours, we can express the Casimir energy by the 
following equations. 

£Casimir(» = £ S ea(^) + Abound- (%) = ~ (s a ky(x) + £bound+ (x)) (21) 
= ~(e sea (x) +£boimd-(z)) _ ^( £ sky(x) + £bound+(^)) (22) 
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e« fi(-v) 




'', -0.2 - / \-02- / 

\ , •' \ , * 

' -0.4 - ; -0.4 - 

-0.6 t ' ' -0.6 - 

Figure 5: The energy densities as a function of x for /i = 4m at 8a = n. The left graph: the solid (dotted) 
line shows the energy density of the negative (positive) continuum states and dashed (dotdashed) line shows 
the energy density of the negative (positive) bound state. The right graph: the solid line shows the sum of 
the negative bound and continuum energy states and the dashed line shows the sum of the positive bound and 
continuum energy states. The Casimir energy density can be represented by the solid line. 



Now we can consider the effect of the Casimir energy on the total energy of the system with a real 
fermion present at the ground state level. For such a system the total energy is the sum of the soliton 
mass and the energy of the single fermion, if the Casimir energy is not taken into account. Since the 
ground state of the system is distorted by the presence of the soliton, this portion should be added. 
Note that for those values of \x for which the energy of the fermion is negative, we should not add this 
energy, since this portion has already been taken into account by the Casimir energy. Figure (JB]) shows 
the total energy (solid line) along with the Casimir energy (dashed line), the soliton mass (dotdashed 
line) and the energy of the fermion in the ground state (dotted line). In this figure we just depict 
the binding energy as long as it is positive. The main portion in the total energy originates from the 
soliton mass. Both the fermion and Casimir energies are multiplied by 10 to make them more visible 
on the graph. 




Figure 6: The total energy (the sum of the energy of a real fermion, the soliton mass and the Casimir energy) 
as a function of [i. The solid line shows the total energy, dashed line the Casimir energy, dotdashed line the 
energy of the soliton [M c \ given in Eq. (J3J)), and the dotted line the energy of a real fermion. The fermion 
energy and Casimir energy are multiplied by 10 to make them more visible. 
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4 Conclusion 



In this paper we have investigated a system consisting of a pseudoscalar field interacting non-linearly 
with a Fermi field. The pseudoscalar field in this model is static and prescribed, and by varying its 
parameters, that is the slope \x and its asymptotic values ±60, we can form different configurations 
including the ones with a proper topological charge. Using the exact fermionic spectrum for this 
model, we have calculated the exact value of the Casimir energy by subtracting the vacuum energy 
in the presence and absence of the soliton. We have displayed the Casimir energy as a function of 
#0 and [i and discovered that this energy has a sharp maximum at the points where the Fermi field 
has zero-mode bound states which would correspond to a fractional fermion number for the soliton if 
particle conjugation symmetry was not broken. Therefore, our analysis shows that the configurations 
which have zero fermionic modes are less likely to occur. As expected, the Casimir energy tends to 
zero as \x — > 00 when 9 is an integer multiple of ir and as fi — > for all values of 60 ■ Furthermore, we 
have shown the Casimir energy as a function of 60 at fi = Am. This graph has two mild cusps at the 
points where a bound state crosses E = 0, and this is the generic behavior of the Casimir energy for 
any values of [i. Then, we have exhibited the contributions of the bound and continuum states to the 
total energy densities for the states with E < and E > 0, separately and found out that they are 
the mirror images of each other at each point of space. Finally, we have added the Casimir energy to 
the total energy of a system consisting of a soliton and a real fermion in the ground state. We have 
concluded that although the contribution of the Casimir energy is small compared to the soliton mass, 
it is not negligible in comparison with the fermion energy. 
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